Eigenvector statistics in non-Hermitian random matrix ensembles 



00 ■ 

a^ ■ 
o^ : 

Oh! 

(D . 
00 ■ 



J. T. Chalker and B. Mehlig 
Theoretical Physics, University of Oxford, 1 Keble Road, 0X1 3NP, United Kingdom 

(February 1, 2008) 

We study statistical properties of the eigenvectors of non-Hermitian random matrices, concen- 
trating on Ginibre's complex Gaussian ensemble, in which the real and imaginary parts of each 
element of an A'' x A'^ matrix, J, are independent random variables. Calculating ensemble averages 
based on the quantity (Lc.|L^) (RfjlRa) , where {Lc,\ and |_R;3) are left and right eigenvectors of J, 
we show for large A'^ that eigenvectors associated with a pair of eigenvalues are highly correlated if 
the two eigenvalues lie close in the complex plane. We examine consequences of these correlations 
that are likely to be important in physical applications. 
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An understanding of statistical properties of ensembles 
of random matrices has proved useful in many different 
areas of physics . Because the first application of ran- 
dom matrix theory, and one still of great importance, was 
to represent the Hamiltonian of a non-integrable quan- 
tum system, the early work of Wigner, Dyson and oth- 
ers focussed on ensembles of real symmetric or complex 
Hermitian matrices. Eigenvector distributions in these 
ensembles are of limited interest, being determined by 
the Haar measure on the group that leaves the ensemble 
invariant. Instead, the concern is mainly with eigenvalue 
correlations, about which a great deal is now known [Q. 

More recently, the spectral properties of random non- 
Hermitian operators have attracted attention in a variety 
of contexts, including: neural network dynamics the 
quantum mechanics of open systems |^]; the statistical 
mechanics of flux lines in superconductors with colum- 
nar disorder ; classical diffusion in random media ; 
and biological growth problems px[ |. The corresponding 
ensembles of real asymmetric and general complex ma- 
trices were first studied by Ginibre , Girko [|l2| , and 
Sommers and co-workers [|[^ . The eigenvalues in these 
ensembles are, of course, not restricted to the real axis, 
but rather distributed over an area in the complex plane. 
Their density and correlations have been investigated in 
considerable detail 

By contrast, eigenvector statistics in non-Hermitian 
random matrix ensembles have not, so far as we know, 
previously been examined. The existence of distinct sets 
of left and right eigenvectors means that invariance of the 
ensemble, under 0{N) or U{N) transformations as ap- 
propriate, is a rather weak constraint on the joint eigen- 
vector distribution: it generates no information on the 
relative orientations of the two sets of vectors. We show 
in this paper that there are, in fact, remarkable correla- 
tions between left and right eigenvectors. These correla- 
tions are likely to be important in physical applications 
of non-Hermitian random matrix ensembles. We illus- 
trate their significance by discussing two consequences, 
involving extreme sensitivity of spectra to perturbations, 
and transients in time-evolution, which are recognised in 
other contexts as typical of non-normal operators p5|-p7[ . 



We consider, following Ginibre [Q, the Gaussian en- 
semble of general complex N x N matrices, J, having 
independent matrix elements, Jfcj, distributed with prob- 
ability 



N 

P{J) dJ (X exp(-iVTr [ J jt] ) [] dj^i dJ'^i 

fe,;=i 
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where Jki 



-f iJ'' with J' 



kl 



and J'^i 



real. Denot- 



ing ensemble averages by (• • •) and complex conjuga- 
tion with an overbar, the only non-zero cumulant of J 
is {Jki Jki) = 1/N. 

The eigenvalues, Aq, of J are distributed in the com- 
plex plane with, in the limit N —^ oo, constant density 
inside a disc of unit radius, centered on the origin. They 
are non-degenerate with probability one, and in this case 
the left and right eigenvectors, (L^l and \Ra), which sat- 
isfy 



J \Ra) — Aa \Ra) 
{La I J — {La I Aa 



(2) 



form two complete, biorthogonal sets, and can be nor- 
malised so that 



{LalRp) 



(3) 



We indicate Hermitian conjugates of vectors in the usual 
way, so that, for example, \La) satisfies \La) = 

Act \La)- 

We investigate eigenvector correlations mainly by cal- 
culating ensemble averages of combinations of scalar 
products. Noting that Eqs. (|^) and are invariant un- 
der a scale transformation \Ra) Ca \Ra) and (L/jj 
(L/jI C^^, one recognises that only those combinations in- 
variant under this transformation should be considered. 
The simplest such combination, involving two eigenvec- 
tors, is fixed by Eq. (^) ; the simplest non-trivial quantity 
is thus the matrix of overlaps 



Oaf! — {La\Lp) {RplRa) 



(4) 
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and we shall focus on this throughout the paper. It is 
convenient to define local averages of diagonal and off- 
diagonal elements of the overlap matrix, 



0{z) 



\iV ■ 
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Correspondingly, the density of states is defined as d{z) — 
{^~^ Ea ^(-^ - ^"))- In ttie limit N oo, d{z) = tx-^ 
for 1^1 < 1 and d[z) — otherwise 11|. 

We have been able to obtain exact expressions for 
0{zx) and 0(zi,Z2). For iV > 1, |zi - Z2I 7^ and 
jzil, |z2| < 1 these simplify to 



0(zi)^^(l-|zi|^ 



0(2:1, Z2) 



1 1 - Z1Z2 

'^W-Z2\^ 



(7) 
(8) 



For |zi|, |z2| > 1, both densities vanish as — > 00. To 
display the form of 0(zi, Z2) as \z\ — Z2 \ ^ 0, it is neces- 
sary to express |zi — Z2I in units of the separation between 
adjacent eigenvalues, introducing z+ ~ (zi -t- Z2)/2 and 
Lo = yfN (zi - Z2)- We obtain, for |z+| < 1, < \fN 
and > 1 



0(Z1,Z2)--7V' 



(l-(l + |c.ne-H^) . (9) 



Eqs. (0) ~ (|9|) constitute our main results. Before outlin- 
ing our derivation, we discuss their significance. 

First, we stress the dramatic difference between the 
behaviour of Oap in this general complex ensemble and 
its behaviour in the case of Hermitian matrices, for which 
Oaf3 — 5af3- The fact that, by contrast, Oaa ^ N in the 
non-Hermitian ensemble can be understood as the be- 
haviour which results if {La\ and \Ra) are independent 
random vectors, subject to the normalisation of Eq. (|^). 
Moreover, large values for the diagonal elements of the 
matrix Oaf3 must be accompanied by some large (or many 
small) off-diagonal elements, since the two are linked by 
a sum rule that follows from completeness. 



Indeed, Eq. (^) implies 

O^f, ^ 0{Z^, Z2)/(^ ^(^1 " ^z^) ^(^2 - A,)) (11) 



and hence, from Eq. (^, Oap ^ —N if Aq and A^ are 
neighbouring eigenvalues in the complex plane, so that 
(typically) lu ^ 1. 



An immediate consequence of large values of Oaa is 
that the spectrum of J for a given realisation has ex- 
treme sensitivity to perturbations. Consider for definite- 
ness J = 008(6*) Ji -I- sin(0) J2, where 9 is real and Ji and 
J2 are both drawn independently from the ensemble of 
Eq. (|^), so that J moves through the ensemble as 9 varies. 
Then 



\dXa/d9f ^\{La\dJ/d9\Ra)\\ 

Performing the ensemble average, and using Eq. 
obtains for the mean square eigenvalue velocity 



{\dX/d9\^) = 



1 



1 



|AP 



(12) 



(13) 



This result should be contrasted in magnitude with the 
analogous one for Hermitian matrices ||l^. Let H = 
cos{9)Hi +sm{9)H2, where Hi and H2 are complex Her- 
mitian N X N matrices, drawn independently from the 
Gaussian unitary ensemble, in which the non-zero cumu- 
lants are {Hu Hik) = {HuiHu) = l/N. Let E be an 
eigenvalue of H. Then for N ^ oo, ~\ < E < \ and 
{[dE/d9Y) = l/N. Thus the eigenvalues of the A^ x A^ 
random non-Hermitian matrix are 0{N) times more sen- 
sitive to perturbations than those of the Hermitian ma- 
trix. Such sensitivity is known to be a typical property of 
non-normal operators |l^]. Despite the sensitivity of in- 
dividual eigenvalues to perturbations, it is reasonable to 
expect some stability in the structure of the spectrum as 
a whole, since the perturbations considered merely take 
a random matrix from one realisation to another. Such 
stability arises from the fact that, although the mean 
square velocity of Eq. ( p^ is large for eigenvalues within 
the unit disc, it vanishes as the boundary to the support 
of the spectrum is approached. Conversely, anticipating 
this stability, we have a rationalisation of the fact that, 
from Eqs. (0),(|) and (|ll]), the amplitudes of the 0{N) 
contributions to Oaa and Oap vanish as |zi| 1. 

The large off-diagonal elements of Oap are significant 
in situations in which J is the generator of evolution in 
real or imaginary time. Settings of this type represent 
one of the main physical applications of non-normal oper- 
ators P,^,p|,p^ . To be specific, consider a model problem 
in which 



(10) so that 



|Kt)) = (j-i)K<)) 



|«(t)) =^|i?«)/((A„)(L„|7/(0)) 



(14) 



(15) 



with /t(A) = exp ( [A — l]i), where we use (J — 1) rather 
than J in Eq. (0^ for convenience, to suppress exponen- 
tial growth. Ensemble averaging with (it(0)|u(0)) = 1 
leads to 



{{u{t)\u{t))) 



a(3 



l/t(A/5) 



(16) 



2 



and for < ^ 1 and iV — > oo we find [fLQl 



{{u{t)\u{t))) ^ (Ant) 



-1/2 



(17) 



This behaviour should be compared with the much faster 
decay that would result from the same spectrum if the 
eigenvectors were orthogonal. In the same regime, the 
replacement Oap — > S^p transforms Eq. (^h into 



-1/2 



(18) 



Thus, eigenvector correlations may be as significant as 
eigenvalue distributions in determining evolution at in- 
termediate times, a fact of established importance in hy- 
drodynamic stability theory |]l5|,^ . 

Finally, it is interesting to ask about, not only the av- 
erage behaviour of the overlap matrix, but also its fluctu- 
ations. In fact, Oap is typically large if the matrix J has 
an eigenvalue which is almost degenerate with or A^, 
and as a result, the probability distribution of Oai3 has 
a power-law tail extending to large |Oq./3|. To illustrate 
this, we consider N = 2, for which we can calculate ex- 
actly the probability distribution, P(Paa), of a diagonal 
element of the overlap matrix. We find 



P{0^c.) = 4 



e(Oc»a- 1) 

(2 0„„-l)3 



(19) 



where Q{x) — 1 for x > and zero otherwise. This im- 
plies in particular that the second and higher moments 
of Oaa diverge. We expect, from Eq. (po|), below, similar 
behaviour for TV > 2 and (provided N > 2) for Oafi with 
a ^(3. 

In the remainder of this paper we sketch our calcula- 
tions and show how the results summarized above can 
be generalized. Calculations for the ensemble of Eq. (|l|) 
can be done by extending the classical methods of Dyson 
and Ginibre, while more general problems are most con- 
veniently treated via ensemble-averaged resolvents, using 
the techniques of Refs. [|l]j2|J|. 

A direct computation of averages of Oaf} involves a 
2A^^-fold integration over the complex matrix elements 
Jki- The integral is simplified considerably by chang- 
ing variables as described in pO[ |. Reducing J by a uni- 
tary transformation, U, to upper triangular form, so that 
T = UlJU has Tki = for k > I, we use as N{N + 1) 
coordinates, the real and imaginary parts of the non-zero 
elements of Tki , and take the remaining coordinates from 
U itself. The required Jacobian is given by Mehta |2^ . 

In this basis, the diagonal elements of T are the eigen- 
values, Tkk = Afc. The first two pairs of eigenvec- 
tors are \Ri) = (1, 0, . . . , 0)t, {Li\ = (1, 62, 63, . . . , ^at), 
\R2) = (-62,l,0,...,0)t and (L2I = (0, 1, dg, • • • , dw), 
where the coefficients bi and di are determined by the 
recursion relations 



Ai 



p-i 

y. bqTqp , 
^ q=l 



Ao - A, 



p-1 



with 61 = 1, di = and d2 = 1. Correspondingly, the 
overlaps are 



N 



1=1 

N 

O12 = -h E ■ 



(21) 
(22) 



Performing the integrals over U and Tki with k < I, 0{z) 
and 0(zi,Z2) can be expressed as averages with respect 
to the joint probability density of the eigenvalues 



N 



P(Ai,..,Aw)«exp(-ivE|Afen Y[\X:-\jf- (23) 

fc=l l<i<j<N 

Defining (• • •)p as an average with ([2^), we find 



2<j<N 



and 



0(21, Z2) = - l)(^d{zi - Ai) 5iz2 - A2) 

n (1 



(25) 



iV|Ai-A2|2 



3<j<JV 



iV(Ai-A,)(A2-A,)//p 



Performing the integrals over eigenvalues in Eqs. ( p4[ ) and 
(psi), we obtain explicit expressions in terms of N x N 
determinants. 

For zi — 0, we are able to evaluate these determinants 
in closed form by recursion, and to simplify the result 
further for iV 3> 1. For zi 7^ we are forced to take a 
less direct approach, which numerical tests show is a good 
approximation for finite N, and which we can prove is ex- 
act in the limit — > 00. We separate the contributions 
to each of the Eqs. ( p^ and ( p5| ) into two factors: one 
from the M eigenvalues closest to zi, with M ^ 1, and 
another from the remaining eigenvalues. The first can 
be evaluated using our result for zi = 0, while the sec- 
ond can be calculated neglecting eigenvalue correlations, 
because its fiuctuations vanish as M — > 00. Their com- 
bination is independent of AI, for M large, as it should 
be, and is as displayed in Eqs. (0) - (^. 

An entirely different approach is necessary in order to 
treat other random matrix ensembles with ease, or to de- 
velop approximation schemes for spatially extended prob- 
lems such as those of Refs. [p|-^,p|To . For these purposes 
we take as central objects the ensemble averages of prod- 
ucts of the resolvents, (zi — J)~^ and (z2 — J'')~^- As a 
demonstration, we examine a matrix ensemble with the 
probability distribution 



P( J) dJ 



(26) 



N 



N 



9=1 



exp ( - Y^Tr [ J Jt - r Re J J] j J] d4i d^i , 



k,l=l 



3 



with r real and — 1 < t < 1 . The non-zero cumulants are 
{Jki Jki) = and {Jki Jik) = t/N. This distribution, 
introduced in ||] , interpolates between the Gaussian uni- 
tary ensemble of Hermitian matrices, for r = 1, Ginibre's 
ensemble (Eq.^) for t — 0, and complex antisymmetric 
matrices for t = —1. In the limit N ^ oo, the eigenvalue 
density has the uniform value d{z) — [7r(l — t^)]^^ within 
the ellipse defined by [Re z/(l - r)]^ -I- [Imz/(1 + t)]2 < 1 
and is zero elsewhere 

We treat the ensemble using the techniques de- 
scribed in 1^,^,^. These generate an expansion for 
0{zi,Z2) in powers of {zi — Z2)/N, and hence give 
0{zi,Z2) exactly in the limit N — > oo, but supply in- 
formation about 0{z) only indirectly, via the sum rule of 
Eq. 10 , We start by considering the 2N x 2N Hermitian 
matrix H = Hq -\- Hi, 



D{zi,Z2) 



(34) 



(1-^2)2, (^^^2)^^^^ ^^(^2^^^) 

(l-r2) 

for zi and Z2 within the support of the density of states, 
and zero otherwise. Since we have taken zi 7^ Z2, this is 
simply 0(zi,Z2), and for t = 0, Eq. (||) is reproduced. 
For 1— r ^ 1, on the other hand, 0(zi, Z2)/[(i(zi)(i(z2)] cx; 
(1 — r), so that 0{zi, z^) /[d{zi)d{z2)\ vanishes in the Her- 
mitian limit r — > 1, as expected. 
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Ha 



, Hi 



z-J 



with real ?7 > 0, and its inverse 

Gil G12 



G = H-^ = 



G21 G22 



(27) 



(28) 



The resolvents are obtained taking 77 ^ 0: in this limit, 
G21 = (z — J)^^ and G12 = (z — J^)^^- Expanding the 
Green function G as a power series in in Hi, its ensemble 
average can be written 



(G) = Go + GoS(G> 



(29) 



where Go — Hq and S is a self-energy. In the limit 
N ^ 00, the self-consistent Born approximation for S 
is exact 22 |], and the eigenvalue density can be ob- 
tained as 



d{z) = lim (^TrG2i(z; 

TT OZ V^O \ N 



(30) 



To study eigenvector correlations, it is necessary to cal- 
culate averages of products of G's. In particular, the 
density 



-D(zi, Z2) = S{zi - Z2)0{zi) -f 0(zi, Z2) 
can be written as 

1 d d 

■K^ dzi dz2 
I 1 



-D(zi,Z2) 



X lim(-TrG2i(zi)Gi2(z2)\ 

r;^0 \ iV / 



We therefore calculate il(zi,Z2) = (G(zi 
which obeys a Bethe-Salpeter equation 
R^{zi, Z2) = (Go(zi)) ® (Go(z2)), 

R = Rq + RqYR . 



(31) 



(32) 



Giz2)), 
writing 



(33) 



In the limit N ^ 00 with zi 7^ Z2, the vertex is simply 
r = diag(l, T, T, 1). Solving Eq. (p3|) for R, we obtain 
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